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Polysilicon fuses are one time programmable memory elements which allow the calibration of integrated circuits at 
wafer and package level. We present a zero dimensional lumped parameter model of the programming of fuses made 
from a combination of tungsten silicide and polycrystalline silicon. The components of the model are an electrical 
model, a thermal model and a flow model. The electrical model describes the temperature and geometry dependent 
resistance of the fuse. The thermal model describes the heating and melting of the fuse and its surroundings. The 
flow model describes the disconnection of the fuse by electromigration driven flow of silica. The model generates 
quantitatively accurate results and reproduces trends with applied voltage and fuse size. 
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Polysilicon fuses are one time programmable memory el- 
ements in silicon chips. They are utilised as permanent on 
chip memory storage for calibration factors in precision ana- 
logue circuits and can be programmed (blown) at both wafer 
and package levels. Polysilicon fuses are constructed from 
layers of polycrystalline silicon (polysilicon) and metal sili- 
cides. Theories of the blowing process which have been put 
forward invoke grain recrystallisation^, melting and vapouri- 
sation^ and electromigration-'. 

One fuse design is shown in Figure [T] The fuse is com- 
pletely solid state and consists of layers of conducting tung- 
sten silicide and semiconducting polysilicon surrounded by 
insulating silica. The thickness of the polysiUcon and tung- 
sten silicide layers are ai and 02 respectively, the width of 
the fuse is h and its length is c. The electrodes on either side 
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FIG. 1. Polysilicon fuse geometry, (a) Top view, (b) Transverse 
cross-section. 
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of the fuse are made from the same materials with the same 
thickness but are much larger in the other two directions. 

The fuse is blown by the application of a fixed potential 
bias across the fuse. Current versus time graphs are shown 
for different voltages and different size fuses in Figure |2] The 
microstructure of a blown 0.3 /^m x 0.3 /im x 1.5 /im fuse 
is shown in Figure [3] The following features of the blowing 
process are apparent: (i) there is evidence of melting of the 
fuse, and of parts of the electrodes and the silica surrounding 
the fuse; (ii) the fuse is disconnected by both a void and by 
a layer of silica. A model of the fuse blowing process must 
explain these observations. 

We propose that the fuse blowing process occurs by the fol- 
lowing steps. (1) When the current starts to flow, it flows 
mostly through the metallic silicide. As the fuse heats up, 
the sihcide conductivity reduces but the conductivity of the 
semiconducting polysilicon increases until it is the dominant 
conductor The overall effect is that the conductivity of the 
fuse first increases then decreases. (2) Eventually the fuse, 
parts of the electrodes and some of the silica surrounding the 
fuse melt. (3) Electromigration drives the fuse material to- 
wards the cathode. This displaces the insulating molten silica 
towards the anode, where it begins to pinch off the fuse from 
the anode. (4) When pinch off is nearly complete the fuse 
ruptures forming the void (this process is driven by the con- 
traction of silicon on melting). 

To model the fuse blowing process we use a lumped param- 
eter approach. This model builds on work reported in Ref. 4l 
The lumped parameter model of the fuse consists of ordinary 
differential equations for the temperature, T, of the fuse, and 
the height, h, and width, w, of the pinched off region. These 
variables are illustrated in Figure 2] An electronic model de- 
scribes the temperature and geometry dependent resistance of 
the fuse; a thermodynamic model describes the heating and 
melting of the fuse and its surroundings; a flow model de- 
scribes the electromigration driven pinch off of the fuse. The 
rupture of the fuse is not modelled explicitly: it is assumed 
that rupture occurs just as pinch off is complete. 

a. Electronics Electrically the system consists of the 
fuse in series with a transistor The fuse blowing process 
is initiated by a decrease of the transistor resistance from a 
very large value to i?t. This decrease occurs on a nanosecond 
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FIG. 3. Electron micrograph showing the microstructure of a blown 
0.3 X 0.3 x 1.5 fuse. 




FIG. 4. Fuse variables: the temperature T of the fuse and the width, 
w, and height h of the pinched off region. 



FIG. 2. Current time curves for: (a) 6 V applied to a 0.3 fim x 
0.6 /im X 3 fim fuse; (b) 5 V applied to a 0.3 fim x 0.6 /im x 3 /im 
fuse; (c) 4.8 V applied to a 0.3 fim x 0.3 /im x 1.5 /im fuse. 
Grey lines, experimental data. Solid black line, simulations includ- 
ing superheating. Dashed black line, simulations without superheat- 
ing. Experimental curves in (a) and (b) were calculated from data in 
Ref.|2|. 



timescale and is assumed to be instantaneous in the model. 
The fuse can be modelled as two resistors in series, with a to- 
tal resistance Rf. The first resistor represents the part of the 
fuse unaffected by flow and has a temperature dependent re- 
sistance. The second resistor represents the region of the 
fuse in which pinch off occurs and thus has a temperature and 
geometry dependent resistance, i?2- 
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The current / flowing through the fuse is calculated from 
Ohm's law. 

Following Ref. 5, conductivity, (j{T), is taken to be a 
quadratic function of temperature. There is insufficient data 



to separate out the contributions of the polysilicon and tung- 
sten silicide layers, therefore a single conductivity function is 
used for both layers of the fuse. 



(7{T) = (To 1 - a (T - Tn, 
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As the fuse heats up, its conductivity increases until 
reached and then it decreases. 

b. Thermodynamics The temperature, T, of the fuse is 
assumed to be homogeneous throughout the fuse. The equa- 
tion modelling its evolution is 



vC{T) 



dT 
dF 



k3s{T-To) 



(4) 



The left hand side describes the heat capacity of the fuse. The 
terms on the right hand side describe Joule heating and Newto- 
nian cooling respectively. C{T) is a the volume averaged heat 
capacity per unit volume of the fuse calculated from: Ci and 
C2, the heat capacity per unit volume of polysiUcon and tung- 
sten silicide (assumed to be the same in the solid and liquid); 
and AHi and AH2, the enthalpies of fusion per unit volume 
of polysilicon and tungsten silicide. In solving the equations, 
the melting phase transitions (at Tmi and Tm2 in polysilicon 
and tungsten silicide respectively) are assumed to have a finite 
width of ST. V is the volume of the fuse, is the thermal con- 
ductivity of silica, Tq is room temperature and s is the surface 
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area of the fuse, 2c (ai + a2 + b), d is a thermal conduction 
length. A better fit to the data is obtained by allowing su- 
perheating of the tungsten silicide before it melts: this seems 
physically reasonable given the small size of the fuse and the 
rapid heating rates. Superheating is included within the model 
by increasing Tin2 from its equilibrium value of 2300 K. 

c. Fluid mechanics The body force driving electromi- 
gration is given by / = gUqlRf/c where Ug is the charge 
density of charge carriers and g is a constant of proportional- 
ity. The force acts on the conducting silicon and silicide and is 
directed from the anode to the cathode. The viscosities of the 
molten silicon and silicide are much lower than that of molten 
silica. Therefore flow is limited by the motion of silica. The 
electromigration body force sets up a pressure gradient within 
the fuse which is transmitted to the molten silica surrounding 
the fuse and causes it to move with a velocity u = fx^/\2rj 
where -q is the dynamic viscosity of molten silica and x is the 
width of the molten silica region surrounding the fuse. 

The equations describing the pinch-off of the fuse are 

Ah ux fh'^ 

757 = (5) 

dr w rjw 

dw fh^ 
dr rjw 

The first equation expresses conservation of volume, the sec- 
ond is a crude balance between work done by body forces and 
work done against friction in the silica. Flow is only allowed 
to occur once Tn,2 has been exceeded. To allow the equations 
to be solved numerically h and w are initialised with small 
non-zero values. 

The equations were solved numerically using a 4th order 
Runge-Kutta scheme with adaptive time-stepping. Parameter 
values used in the simulations are summarised in Table J] The 
results are shown in Figure |2] for models using the equilibrium 
melting value of Tm2 and the values including superheating 
given in Table |T] Considering the simplicity of the model the 
results are remarkably good and trends with voltage and fuse 
size are correctly predicted. 

A zero dimensional, lumped parameter model of the pro- 
gramming of polysilicon fuses successfully reproduces the 
available experimental data for one type of fuse. The model 
captures trends with voltage and size. This strongly suggests 
that the assumptions used to construct the model are funda- 
mentally correct. To take the model further, more experimen- 
tal data are required, describing the temperature dependent 
conductivity of the individual polysilicon and tungsten sili- 
cide layers for example. This will allow three dimensional 
modelling of the partial differential equations describing the 
system. 
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TABLE I. Values of parameters used in simulations. 



Parameter 


Value 




R, 


5 n 






1700 K 




V 


10 Pas 




riq 


3.16 X 10"' C m" 


-3 


9 


0.2 




To 


3.33 n-^ 


m-1 


Tniax 


1400 K 




To 


300 K 




a 


6.89 X lO"'' K"^ 




Ci 


1.66 X 10^ Jm~ 


3^-1 


Ca 


3.09 X 10® J m- 


3^-1 


AHi 


4.15 X 10^ J m" 


3 


AH2 


8.30 X 10^° J m" 


3 


ST 


50 K 




k3 


1.5 Wm 


-1 g-1 



Property 


(a) 


(b) 


(c) 


0.1 


0.15 /im 


0.15 nm 


0.15 /im 




0.15 fim 


0.15 ^m 


0.15 /im 


b 


0.6 fim 


0.6 ^m 


0.3 /im 


c 


3 fim 


3 /im 


1.5 /im 


d 


0.3 fim 


0.3 /im 


0.15 /im 


X 


0.1 fim 


0.1 /im 


0.08 /im 


V 


5 V 


6 V 


4.8 V 


Tm2 


2400 K 


2500 K 


2500 K 
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